Cold-atoms on a two-dimensional square optical lattice with an alternating potential 
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The cold-atom on a two-dimensional square optical lattice is studied within the hard-core boson 
Hubbard model with an alternating potential. In terms of the quantum Monte Carlo method, it is 
shown explicitly that a supersolid phase emerges due to the presence of the alternating potential. 
For the weak alternating potential, the supersolid state appears for the whole range of hard-core 
boson densities except the half-filling case, where the system is a Mott insulator. However, for the 
strong alternating potential, besides the supersolid and Mott insulating states, a charge density 
wave phase appears. 
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The optical lattice, which is produced by using laser 
beams, allows one to tune relevant parameters (such as 
the confining potential, particle density, interaction, lat- 
tice geometry, lattice constant, etc.), and therefore to 
make it an ideal test-ground for concepts and theories on 
condensed matter physics, especially strongly correlated 
many-body quantum lattice models [U, [3i H I3j [E| ■ Theo- 
retically, it has been argued that the physical properties 
of bosons trapped in an optical lattice can be described 
within the boson Hubbard model. This boson Hubbard 
model and its extended version have been studied ex- 
tensively [i, 0, i, S 113, im, d, and the results show a 
number of interesting features, such as the existence of 
superfluid state, charge density wave (CDW) state, and 
Mott insulating state. 

Over fifty years ago, Penrose and Onsager [13] sug- 
gested a novel state with simultaneous solid and super- 
fluid phases, which intrigued in condensed matter physics 
field. We call such state the supersolid state. Re- 
cently, a possible supersolid state in the solid ''He has 
renewed interests in this new matter state p^ . How- 
ever the cold-atom in an optical lattice is a more promis- 
ing realization of supersolid phase. Theoretical studies 
[1, 0, B S [13, im, have shown that next-nearest- 
neighbor interaction, frustration of lattice, and soften- 
ing of the on-site interaction are favorable for the real- 
ization of the supersolid state. However, based on the 
two-dimensional (2D) hard-core boson Hubbard model 
on a square lattice with nearest-neighbor repulsion in- 
teraction, it has been found that there is no supersolid 
phase due to the soHd-supcrfluid phase separation 0. 
Motivated by the ionic Hubbard model [l3|, whether or 
not the phase diagram of the 2D hard-core boson Hub- 
bard model can be infiuenced by an alternating potential 
is an exciting issue. 

In this paper, we study this important problem. 
Within the hard-core boson Hubbard model on a square 
lattice with an alternating potential, we employ the quan- 
tum Monte Carlo (QMC) method, and show explicitly 
that the supersolid phase emerges due to the presence 
of the alternating potential. In particular, our results 
also indicate that even very weak alternating potential 
can prevent the occurrence of the solid-superfluid phase 



separation, and then induce the supersolid state. 

The 2D hard-core boson Hubbard model on a square 
lattice with an alternating potential can be expressed as 

H^-tY, (alaj + a]ai) - ^Y^n, + A ^(-l)^-+*«n„(l) 

<ij> i i 

where < ij > indicates the sum over nearest neighbor 
sites, the operator a| (a^) creates (destroys) a hard-core 
boson on site i, hi — a|ai is the hard-core boson num- 
ber operator, fi is the chemical potential and therefore 
controls the filling number in the trap, hopping t de- 
scribes the coherent hopping between nearest-neighbor 
sites, and (— 1)(*=:+*!') A is the alternating potential. In 
the hard-core limit, the Hubbard constant [/ — s- oo. In 
this case, each lattice site can be occupied by or 1 hard- 
core boson. Since hard-core bosons are restricted in this 
Hilbert subspace, operators a| and obey commutation 
relations [ai,a]] = at sites i ^ j, while they satisfy 

anti-commutation relations {a^, aj} = 1 on sole site i. 

In the following discussions, our QMC simulation is 
based on the stochastic series expansion method with di- 
rected loop updates [l^]. In Hamiltonian (1), the dif- 
ferent phases are characterized by static staggered struc- 
ture factor 5(Q) with wave vector Q = [tTjTt] and su- 
perfluid density ps, where the structure factor measures 
the diagonal long-range order of the system, and the su- 
perfluid density measures the off-diagonal long-range or- 
der. We emphasize that in the square lattice magnetic 
system, the wave vector Q is an antiferromagnetic wave 
vector. However, in the present square lattice hard-core 
boson system, this wave vector Q is similar to that in 
the magnetic system, and therefore the structure factor 
5(Q) measures the diagonal long-range order with the 
checkerboard pattern. In the present case, the structure 
factor, the superfluid densit y, a nd the hard-core boson 
density are expressed as [l2l.1l7j: 
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where W is the winding number of the bosonic world 
hnes. Due to the presence of the alternating potential 
(_1)(«x+s)a in Hamiltonian (1), there are two sublat- 
tices A and B in the system. This leads to the for- 
mation of two energy bands, where the gap is opened 
by the alternating potential. At half-filling (correspond- 
ing to the hard-core boson density p — 1), both energy 
bands are completely occupied by hard-core bosons, so 
all hard-core bosons are localized, with 5'(Q) = and 
Ps — 0, and the corresponding state of the system is 
a Mott insulating state. For the case of the hard-core 
boson density p = 1/2, the sublattice with minus po- 
tential is completely occupied by the hard-core bosons, 
while the other sublattice is empty, and the correspond- 
ing state of the system is CDW state, with S{Q) ^ 
and Ps = 0. The conditions for the supersolid phase are 
that both S{Q) ^ and 7^ simultaneously. The su- 
persolid state is a novel state with simultaneous diagonal 
(CDW solid) and off-diagonal (superfluid) long-range or- 
ders In the following QMC calculations, we discuss 
the structure factor S'(Q) and the superfluid density ps as 
a function of the hard-core boson density p, the hard-core 
boson density p as a function of the chemical potential /z, 
and the ground state phase diagram of the hard-core bo- 
son Hubbard model (1) with an alternating potential in 
the p— A plane, where all the numerical results have been 
performed on a. N — L x L lattices with L = 16, and the 
inverse temperature was chosen as /? = 2L, which is low 
enough for the numerical simulation of the ground state 
properties of the cold-atom in an optical lattice. The 
finite size extrapolation of the following measurements 
converges to a finite value, indicating that the obtained 
phase is stable in the thermodynamic limit. 

Firstly, we plot the superfluid density ps (circle) and 
the structure factor 5(Q) (triangle) as a function of the 
hard-core boson density p for parameter (a) A/t = 1, (b) 
A/t = 2, (c) A/t = 4, and (d) A/t = 6 in Fig. 1, where 
the superfluid density ps and the structure factor S{Q) 
are changed dramatically with the alternating potential 
A. In particular, for any alternating potential, the su- 
perfluid density and the structure factor are symmetrical 
around p = 1/2, reflecting the particle-hole symmetry. 
Moreover, for any alternating potential, the Mott insu- 
lating state appears only at half- filling {p — 1). Surpris- 
ingly, for the weak alternating potential, both S{Q) 7^ 
and Ps 7^ at < p < 1. In this case, CDW state is 
absent, and the supersolid state appears for all hard-core 
boson densities except p = and p = 1 . This is much dif- 
ferent from the case of the 2D hard-core boson Hubbard 
model on a square lattice with nearest-neighbor repulsion 
interaction, where the supersolid phase is thermodynam- 
ically unstable, and therefore no supersolid phase can be 
observed due to the CDW solid-superfluid phase separa- 
tion [l,[l3 • However, for the strong alternating potential, 
besides the Mott insulating and supersolid states, CDW 
state appears at p = 1/2. Moreover, the superfluid den- 
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FIG. 1: The superfluid density ps (circle) and structure factor 
^(Q) (triangle) as a function of the hard-core boson density 
p for (a) A/t = 1, (b) A/t = 2, (c) A/t = 4, and (d) A/t = 6 
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FIG. 2: The hard-core boson density p as a function of the 
chemical potential p for A/t = 1 (square), A/t — 2 (circle), 
A/t = 4 (triangle down), and A/t — 6 (triangle up). 



sity Ps ^ at p = 1/2 while S'(Q) ^ 1/4, which is an 
exact value of the static structure factor of the checker- 
board solid. 

For a further confirmation of the above obtained con- 
clusion, we have made a series of calculations for the 
hard-core boson density p as a function of the chemical 
potential p, and the results for A/t = 1 (square), A/t = 2 
(circle), A/t — 4 (triangle down), and A/t = 6 (triangle 
up) are plotted in Fig. 2. Obviously, for any alternat- 
ing potential, p varies continuously with p., reflecting a 
phase transition of second order. When p is less than 
—A, hard-core bosons begin to reside in the optical lat- 
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FIG. 3: The ground state phase diagram in the /i — A plane. 
Inset: the corresponding result for t = 0. 
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FIG. 4: The compressibility k as a function of the hard-core 
boson density p for (a) A/f = 1, (b) A/t = 2, (c) A/f = 4, 
and (d) A/f = 6, where the unit of k, is 1/t. 



tice due to the presence of hopping t. This critical value 
[ic obtained from QMC simulations is quantitatively con- 
sistent with [ic = — \/A^ -f 16f^ obtained from analytical 
calculations. Moreover, a plateau at p = 1/2 appears for 
the strong alternating potential. In comparison with the 
results in Fig. 1, we therefore find that the state in this 
plateau corresponds to a CDW state. This plateau de- 
creases with the decrease of A, and vanishes in the weak 
alternating potential. This result could be anticipated, 
since the strong alternating potential favors CDW state. 
For the large values of /i, the system is at half-filling, so 
this leads to the formation of an uniform Mott insulat- 
ing state. Furthermore, the range of the supersolid state 
increases with the decrease of A, and then in the weak 
alternating potential, the ground state for all hard-core 
boson densities is the supersolid state except the Mott 
insulating state ai p = 1. This is a remarkable result, 
since it indicates that even very weak alternating po- 
tential can prevent the occurrence of the solid-superfluid 
phase separation, and then induce the supersolid state. 
To show these points clearly, we plot the ground state 
phase diagram in the /i — A plane in Fig. 3. In the 
case of hopping t ^ 0, the hard-core boson Hubbard 
model (1) with an alternating potential is reduced to 
H = —p'Y^^rii + A ^^(— l)*^+'!'ni. For comparison, the 
phase diagram in this case {t = 0) is also plotted in Fig. 
3 (inset). It appears that the supersolid state, which ap- 
peares for t is not observed for t = Q, and the ground 
state is CDW state at p = 1/2 or Mott insulating state 
at p = 1. On the other hand, in the case of the alter- 
nating potential A = 0, the model Hamiltonian (1) is 
reduced to the simple hard-core boson Hubbard model 
H = -t^^ijy{a\aj + a]ai) ~ fJ-J2i "ii where the super- 
fluid phase [S{Q) = and ps ^ 0] appears at < p < 1 
[T^ . and no supersolid state is observed. In this case. 



the supersolid phase in Fig. 3 does not include the line 
of A = 0. These results therefore confirm that for finite 
hopping t, the supersolid state emerges due to the pres- 
ence of the alternating potential in the hard-core boson 
Hubbard model, in particular, for the weak alternating 
potential < A/i < 2, the ground state for all hard-core 
boson densities is the supersolid state except the Mott 
insulating state at p = 1 and the empty state at p = 0. 

One of the characteristic features for the thermody- 
namically stable state is that it has a positive compress- 
ibility. To better perceive the stable states shown in Fig. 

3, we turn now to the discussion of the compressibility of 
the system, which is defined by k = /3iV(p^ — (p)^) [ISj . 
We have performed a QMC simulation for the compress- 
ibility K, and the results of function of p for (a) 
A/t = 1, (b) A/t = 2, (c) A/t = 4, and (d) A/t = 6 
are plotted in Fig. 4. In the case of hopping t = 0, there 
are only three states as shown in Fig. 3: the Mott state 
at p = 1, CDW state at p = 1/2, and the empty state at 
p — 0, then zero compressibility for all these three states 
is found. However, in the present results shown in Fig. 

4, we find that for the strong alternating potential A, 
zero compressibility is observed at p = 1/2, which show 
that CDW phase is an incompressible state. Moreover, 
for all other densities except p = and p = 1, positive 
compressibility is observed, which indicates the absence 
of the phase separation, and therefore confirms that the 
supersolid state is the ground state of the system in this 
density range. 

Now we give some physical interpretation to the above 
obtained results. In analogy to the spinful Hubbard 
model [l^, the hard-core boson Hubbard model (1) with 
an alternating potential describes a competition between 
the kinetic energy (pt) and alternating potential (A). In 
the case of p = 0, the system is empty. With the increase 
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of the number of hard-core bosons (0 < p < 1), there is a 
gain of kinetic energy per hard-core boson proportional 
to t due to the hopping. At the same time, CDW or- 
der is reduced, costing an energy of approximately A per 
site. The kinetic energy pt favors the superfluid state and 
tends to reduce CDW order, while the alternating poten- 
tial A favors CDW order and results in frustration of the 
kinetic energy, then this competition results in the emer- 
gence of the supersolid phase. For the weak alternating 
potential, the kinetic energy of the system is larger than 
the potential energy, therefore CDW order is reduced of 
several orders and the system is supersolid. With the in- 
crease of potential energy, CDW order is enhanced, and 
in the region where the potential energy of the system is 
larger than the kinetic energy, CDW state appears at the 
hard-core boson density p = 1/2. In the case of p = 1, 
both energy bands are completely occupied by hard-core 
bosons, and therefore all hard-core bosons are localized. 
In this case, there is no "charge degree of freedom" in the 
system, which is a Mott insulator. 

In conclusion, within the hard-core boson Hubbard 
model with an alternating potential, we have employed 
QMC method to investigate the ground-state physical 
properties of the cold-atom on a optical lattice. Our re- 
sults show explicitly that the supersolid phase emerges 
due to the presence of the alternating potential. At half- 
filling p = I, the system is a Mott insulator. For the weak 



alternating potential, the supersolid state appears for the 
whole range of hard-core boson densities < p < 1. How- 
ever, for the strong alternating potential, the supersolid 
state can be observed for < p < 1 except p = 1/2, 
where CDW state appears. 

Experimentally, the 2D optical lattice has been real- 
ized by using four beams of lasers [13]. Moreover, the 
cold-atom on this 2D optical lattice has been studied ex- 
perimentally. In particular, a superlattice potential that 
is similar to the alternating potential discussed in the 
present case has been realized experimentally by using 
multiple wavelength laser beams pll [2^. Therefore, it 
is possible that our present results of the cold-atom on a 
two-dimensional square optical lattice with an alternat- 
ing potential can be verified by further experiments. 
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